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The observables ck and AMk play a prominent role in particle physics due to their sensitivity 
to new physics at short distances. To take advantage of this potential, a firm theoretical prediction 
of the standard-model background is essential. The charm-quark contribution is a major source of 
theoretical uncertainty. We address this issue by performing a next-to-next-to-leading-order QCD 
analysis of the charm-quark contribution r\ cc to the effective |AS| = 2 Hamiltonian in the standard 
model. We find a large positive shift of 36%, leading to r\ cc = 1.87(76). This result might cast 
doubt on the validity of the perturbative expansion; we discuss possible solutions. Finally, we 
give an updated value of the standard-model prediction for |ejc| = 1.81(28) x 10 -3 and AAffp = 
3.1(1.2) x l(T 15 GeV. 



Strangeness-changing neutral-current transitions play 
an important role in particle physics. The parameter ck, 
measuring indirect CP violation in the neutral Kaon sys- 
tem, has received increased attention recently due to the 
discrepancy between the theoretical prediction and the 
experimental measurement [1-4 . In addition, together 
with the Kaon mass difference AMk, it provides strong 
constraints on many models of new physics. 

Theoretical predictions for AMfP and €k are calcu- 
lated in the framework of effective field theories, which 
allow to separate short- and long-distance contributions, 
and to sum all terms which are enhanced by powers of 
large logarithms \og{m 2 c /My V ) using the renormalisation 
group. The relevant \AS\ — 2 effective Hamiltonian in 
the three-quark theory reads 

n f=3 2 = ^ M w[^ cc S{x c ) + \ 2 tm S(x t ) 

+ 2\ c \ t ri c tS(x c ,Xt)]b(n)Qs2 + h.c. (1) 

where Gf is the Fermi constant, A^ = V* d Vi S comprises 
the Cabibbo-Kobayashi-Maskawa matrix elements, and 
Qs2 — (sLJud-L) 2 is the leading local four-quark operator 
that induces the IAS*! = 2 transition, defined in terms of 
the left-handed s- and d-quark fields. The parameter b(p) 
is factored out such that 
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where fx is the Kaon decay constant, is a 
renormalisation-group invariant quantity compris- 
ing the hadronic matrix element. It can be calculated 
on the lattice with high precision [5HS]. 

The loop functions S can be found, for instance, 
in [TD] . The QCD and logarithmic corrections are con- 
tained in the 77 factors and are known at next-to-leading 
order (NLO) for the dominant top-quark contribution 
(rjtt = 0.5765(65) [H]). The relative suppression of the 



top-quark contribution by the small imaginary part of 
A 2 , relevant for ck, lets the charm-quark contributions 
compete in size. We have already performed a next-to- 
next-to-leading-order (NNLO) calculation of the charm- 
top contribution (r] ct = 0.496(47) [H]). Here, we focus 
on the charm-quark contribution, known until now at 
NLO, with a substantial error {r] cc = 1.40(35) [31 H5]). 
It multiplies S(x c ) = x c + 0(x 2 ), where x c = m^/M^ 
and m c = m c (m c ) is the MS charm-quark mass. The 
Glashow-Iliopoulos-Maiani (GIM) mechanism cancels a 
potential large logarithm at leading order (LO). 

The charm-quark contribution r\ cc determines the 
short-distance part of the Kaon mass difference AMfP 
and enters e% with a negative sign. The large remain- 
ing scale uncertainty at NLO hints at potentially sizeable 
NNLO corrections; we confirm this expectation by an ex- 
plicit calculation in this Letter. 

Our calculation proceeds in three steps: determination 
of the initial conditions of the Wilson coefficients at the 
electroweak scale, renormalisation-group evolution to the 
charm-quark scale, and matching onto the effective three- 
quark theory. The new result is the three-loop matching 
condition at the charm-quark scale. 

The effective Hamiltonian in the five- and four-flavour 
theory relevant for rj cc reads 
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Here the current- current operators are given by Q q ± = 

{{ni^l)®{ii^dp±{ni^i)®{ni^dl))/2, where 
a and /3 are colour indices, and we define the evanescent 
operators in such a way that the anomalous dimension 
matrix is diagonal through NNLO QH E]. The GIM 
mechanism cancels a contribution of the |AS| = 2 Hamil- 
tonian above the charm-quark scale; we verified explicitly 
that mixing of dimension-six into dimension-eight opera- 
tors proportional to A 2 does not occur above the charm- 
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FIG. 1. Sample one-, two-, and three-loop diagrams con- 
tributing to the matching at the charm-quark scale. Loopy 
lines are gluons, and straight lines are quarks. The combina- 
tion c — u arises from the GIM mechanism; q denotes any of 
the quarks u, d, s. 



quark scale. 

We take the initial conditions for C±, obtained by a 
NNLO matching calculation at the electroweak scale, 
from Ref. [T3]. The dimension-eight Wilson coeffi- 
cient does not receive a contribution at the electroweak 



scale [TS]. The running of C± to the charm-quark scale 
can be taken up to NNLO from [14]. 

At the scale fi c — 0(m c ) the charm quark is removed 
from the theory as a dynamical degree of freedom. Re- 
quiring the equality of the Green's functions in both the- 
ories at \x c leads to the matching condition 
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which we use to determine the Wilson coefficient C 
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delined implicitly in ( 10 1 below. Here, angle brack- 



ets denote operator matrix elements between s- and d- 
quark external states. Writing (Qsz) — ?"S2(Qs2)^ and 
(QiQj) = Hj/ / (8TT 2 )dij(Qs2)^°K and expanding all quan- 
tities in powers of a s /{Ait), we find the following con- 
tributions to the matching (a sum over i,j = +, — is 
implied) : 
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The strong coupling constant a s is defined in the three- 
quark theory throughout this Letter, and superscripts in 
brackets denote the order of the expansion in a s . Fur- 
thermore, we expand the charm-quark mass defined at 
the scale fi c , viz. m c (fi c ), about m c (m c ), as in Ref. |14j . 

In order to evaluate the Eqs. ([5]), we compute the fi- 
nite parts of one-, two-, and three-loop Feynman dia- 
grams of the type shown in Fig. [TJ the evanescent op- 
erators in the |A£| — 2 sector have been chosen as 
in [T2]. Our NLO result confirms the calculation by 
Herrlich and Nierste [16] for the first time. The three- 
loop matching calculation yields (we use the notation 
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dyVgj; note also 
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where we defined l c — log(fi 2 /m 2 (fi c )), £ n denotes Rie- 
mann's zeta function of n, and the remaining constants 
are defined in [T7]. This result is new. 

Since the calculation of the NNLO contributions to r) cc 
is quite complex, we checked our results in several ways. 
First of all the calculation of the O(10 000) Feynman di- 
agrams, the renormalisation, and the matching calcula- 
tion, has been performed independently by the two of us, 
using a completely different set of computer programs, 
leading to identical results. On the one hand we use 
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qgraf [TS] for generating the diagrams; the evaluation of 
the integrals is then performed using the program pack- 
ages q2e/ exp / M ATAD pU [19] [20] . On the other hand, 
all calculations have been performed using an indepen- 
dent setup, based on Feynarts ,21], Mathematica, and 
FIRE [22]. 

As a further check of our calculation, we verified that 
the matrix elements are finite and independent of the 
gauge-fixing parameter £. We have also checked analyti- 
cally that rj cc is independent of the matching scales \iw 
/ib, and n c to the considered order of the strong coupling 
constant, by expanding the full solution of the renormal- 
isation group equations about the respective matching 
scale. 

The effective Hamiltonian valid below the charm-quark 
threshold contains only the single operator Qs2- The 
renormalisation-group evolution of the Wilson coefficient 
Cg c 2 is described by the evolution matrix corresponding 
to the anomalous dimension of Qs2'- 

C c s c 2 (fi) = U(^fi c )C c s c M. (9) 

We express the coefficient rj cc in a scale- and scheme- 
independent way as 

Vcc = -^r-,C% 2 ( Mc ) [a s ( Mc )f + K^ 1 (fi c ) . (10) 

TTlc \ m c) 

The remaining scale dependence present in ^ is ab- 
sorbed into 

6 ( M ) = K( M )P + (11) 
where, up to second order in a s , 

KM=1+ j>r-^ + j?(^>)\ (i2) 

and the exponent a + = 2/9 is the so-called magic number 
for the operator Q + (the magic numbers as well as the 
matrices J, comprising the higher-order QCD contribu- 
tions to the renormalisation-group evolution, are defined, 
for instance, in [23]). This scale dependence is cancelled 
by the corresponding scale dependence of the hadronic 
matrix element, order by order in perturbation theory. 
Consequently, our result is independent of \i c up to and 
including terms of 0(a 2 ). 

As a first estimate of the theoretical uncertainty of r\ cc 
we study the residual scale dependence, using three dif- 
ferent methods to evaluate the running strong coupling 
constant |14) . Matching at m c (m c ) and varying /i c be- 
tween 1 and 2 GeV (see Fig. [2]) and \iw between 40 and 
160 GeV we find the following numerical value at NNLO, 

r hc = 1.86±0.53 Mc ±0.07 MW ±0.06 Qs ±0.01 mc , (13) 

where we also display the parametric uncertainties 
stemming from the experimental error on a s {Mz) 



0.1 1M;7: [H] and m c (m c ) = 1.279(13) GeV :25J. The 
dependence on the scale /if, and on nit is completely neg- 
ligible [26]. 
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FIG. 2. a function of fi c , matching at p = m c (m c ) and 

fixing )i w = 80 GeV and = 5 GeV. The LO result is rep- 
resented by the double-dotted line. We also show the NLO 
value of rjcc, with the running a s evaluated either by solv- 
ing the renormalisation group equations numerically (dashed 
line), or by first computing the scale parameter Aqcd, either 
explicitly (dotted line) or iteratively (dash-dotted line) - see 
Ref. [27] for the details. The resulting uncertainty is size- 
able at NLO. The solid lines show the corresponding NNLO 
results; now the ambiguity is almost cancelled, whereas the 
residual scale dependence is still large. 

Varying /i c and \iw in the same range as above, we 
find at NLO 

77 c N c LO = 1.38 ± 0.52 Mc ± 0.07„ w ± 0.02 Qs , (14) 

where the error indicated by the subscript "/i c " includes 
the effect of the three ways of determining a s . We have 
included the parametric uncertainty related to cv s ; the 
error resulting from m c is negligible. 

We find a substantial shift from NLO to NNLO for ?y cc ; 
furthermore, we observe that the NNLO calculation does 
not reduce the residual scale dependence (see Fig. ]2j). 
This reveals a bad convergence behaviour of the expan- 
sion of rj cc in a s , even after having summed all terms pro- 
portional to < log {ml/M^Y and < +1 log {m 2 JM^) n 
using the renormalisation group. Some of the higher- 
order terms leading to the residual scale dependence are 
scheme dependent. In order to show that the large scale 
dependence is not artificial and to gain a better under- 
standing of its origin, we expand the full solution of the 
renormalisation group equation. We find, up to terms 
cubic in a s , 

W(a s (m c )) 2 / 9 = 1 + a s (m c )(0.25 + 0.32L C ) 
+ (a s (m c )) 2 (1.20 + 0.03L,, + 0.22i c + 0.27L 2 C ) , 

where L c = log(m 2 c /M^) = -8.28, L b = log(m^/M^) = 
—5.92, and a s (m c ) — 0.35 at three-loop accuracy. Here 
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we neglect the small terms proportional to L\. This re- 
sult is independent of the renormalisation scale and the 
definition of the evanescent operators and depends im- 
plicitly only on the choice of the renormalisation schemes 
used to determine Bk (cf. Eq. (lOl), a s , and the charm- 
quark mass. The large logarithmic terms proportional 
to powers of L b and L c are summed to all orders by the 
renormalisation group in our full result, but the constant 
term of the NNLO correction is more problematic: it is 
almost twice as large as its NLO counterpart. Such large 
constant parts are expected to lead to a large residual 
scale dependence, as we indeed observe in Fig. [2j 

The convergence of the series can be somewhat im- 
proved by expanding the square of the charm-quark mass 
multiplying r\ cc in Eq. in powers of a s , noting that the 
charm-quark mass receives negative corrections, although 
the effect is not substantial at NNLO. 

As a consequence of the discussion above, we propose 
the following temporary prescription: we take rj cc at /i c = 
m c as the central value, and as the theory uncertainty the 
absolute size of the NNLO correction and the residual 
scale dependence, added in quadrature. This leads to 



'lc 



= 1.87 ±0.76. 



(15) 



Compared to the NLO value ?7^ LO (14), this corresponds 
to a positive shift of approximately 36%. The parametric 
uncertainty is essentially negligible with respect to the 
theoretical uncertainty. 

Finally we study the impact of rj cc at NNLO on 
the prediction of \e K \ and AAfJp. We use the in- 
put values from [53], in particular \V c b\ — 4.06(13) x 
10~ 2 , plus TOt(rot) = 163.7(1.1) GeV [25], m b (m b ) = 
4.163(16) GeV [25], A = 0.2255(7) [29], K e = 
0.923(6) [30], £ s = 1.243(28) [31], r) tt = 0.5765(65) pi], 
B K = 0.737(20) [H][22, r]ct = 0.496(47) PI], in the fol- 
lowing formula (we express fj and p through sin 2/3; for a 
discussion and definitions see [II ITU]): 

\e K \ = K t C e B K \V cb \ 2 \ 2 n[\V cb \ 2 {\ - p)vuS(x t ) 

+ VaS(x c , x t ) - f] cc S(x c )] . (16) 



Using the numerical values given above, we obtain 

\e K \ = (1.81 ± 0.14„ cc ± 0.02„ tt ± 0.07„ ct ± 
+ 0.23 ^.l x 10" 3 . 



05 LD 



± 0-23 paramc tri C ) x 10 3 . (17) 



The first three errors correspond to rj cc , rj t t, r] ct , respec- 
tively. The error indicated by LD originates from the 
long-distance contribution, namely Bk and K e , which 
account for 81% and 19% of the long-distance error, re- 
spectively. Half of the parametric error stems from |V c f,| 
(49%), while all other contributions are well below 20%. 
All errors have been added in quadrature. 

Compared to the prediction using the NLO value , 
\e K LO \ = 1.90(27) x 10~ 3 , this corresponds to a shift of 



approximately —5%, and overcompensates the shift of 
±3% found in [12]. The large perturbative corrections 
are thereby partially mitigated in the observable ck ■ 

Finally we estimate the short-distance contribution to 
AM K . Using [10] 



AM k B = ^f K B K M K M w 
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we find AMfP 



= 3.1(1.2) x 10" 15 GeV, where the central 
value accounts for 89% of the measured value. We ne- 
glected the correction due to top quarks, of the order of 
1%. The error is dominated by r] cc (86%) and Bk (6%). 
Unfortunately, the LD contributions to AMr- are poorly 
known; the discussion in Ref. [33] hints at a positive con- 
tribution. In addition, our calculation shows that also 
the SD contribution cannot be computed as reliably as 
thought previously, and thus the prediction of the total 
Kaon mass difference suffers from large uncertainties. 

We have performed the first NNLO QCD analysis of 
the charm-quark contribution r\ cc to the |A5*| = 2 ef- 
fective Hamiltonian 1-[)^3~ 2 . We confirm the analytical 
results for i] cc obtained at NLO in Ref. [16] for the first 
time. 

The discrepancy between our standard-model pre- 
diction and the precisely measured experimental value 
|e^| cx P = 2.228(11) x 10~ 3 24] could be interpreted as 
a tension within the standard model if we got a bet- 
ter control of the theoretical uncertainty. In view of 
the considerable residual scale dependence and the large 
NNLO shift, sizeable corrections beyond NNLO may be 
expected. 

Given the importance of the observable €k , an effort 
should be made to circumvent these difficulties. We see 
at least two possible ways to proceed: in the short run, 
one could make use of the cancellation of the scheme 
dependence between of the parameter Bk and the effec- 
tive Hamiltonian. One could utilize this scheme depen- 
dence (which would affect the quantities J in Eq. (12)) 
to achieve a better convergence of r\ cc . Recently, new 
lattice renormalisation schemes have been employed in 
the determination of Bk [3 EH] \ they use nonexceptional 
momentum configurations, leading to better control over 
lattice uncertainties. Furthermore, they might lead to a 
better convergence at NNLO, as suggested by the good 
perturbative behaviour of the continuum matching for 
the light-quark masses [3SJ [3B]. We encourage the in- 
vestigation of the effects of these schemes also on the 
convergence of the series for r] cc , in particular, at NNLO. 
In the long run, the possibility of calculating the effects 
of a dynamical charm quark on the lattice might seem 
most promising and should be further studied. 
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